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Abstract: Particle physics models with Natural Supersymmetry are characterized by a 



j-^l superpotential parameter /x ~ nih ~ 125 GeV, while third generation squarks have mass 

~ 0.5 — 1.5 TeV. Gluinos should be lighter than several TeV so as not to destabilize the 

K^ lighter squarks. First and second generation sfermions can be at the tens-of-TeV level 

0> ' which yields a decoupling solution to the SUSY flavor and CP problems. Adopting a 

51^ . top-down approach, we delineate the range of GUT scale SUSY model parameters which 

lO ■ leads to a Natural SUSY mass spectrum. We flnd natural SUSY models to be tightly 

cn ! constrained by the 6 — >• 57 branching fraction measurement while it is also difficult but not 

impossible to accommodate a light Higgs scalar of mass ~ 125 GeV. We present several 

benchmark points which are expandable to slopes and planes. Natural SUSY is difficult 

to see at LHC unless some third generation squarks are very light. The top- and bottom- 

k><( ■ squarks cascade decay mainly to higgsino-like charginos and neutralinos via numerous 

?H ■ posibilities, leading to a rather complex set of signatures. Meanwhile, a linear e'^e~ collider 

- - - operating at ^/s ~ 0.25 — 0.5 TeV would be a higgsino factory and is essentially guaranteed 

a SUSY discovery of the low- lying charged and neutral higgsino states. Since thermal 

neutralino cold dark matter is underproduced, we conjecture that the incorporation of a 

Peccei-Quinn sector or light moduli into the theory will augment higgsino dark matter 

production, possibly together with an admixture of axions. We present rates for direct 

and indirect higgsino dark matter detection for the case where light higgsinos dominate 

the dark matter abundance. 
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1. Introduction 

1.1 Impact of LHC sparticle searches 

The search for weak scale supersymmetry (SUSY) [Q, ^, ^] has begun in earnest at the 
CERN Large Hadron ColHder (LHC). From a non-observation of multi-jet plus multi- 
lepton -|- E^^^^ events with or without tagged 6-jets in a data sample of ~ 1 — 4.5 fb~^, the 
CMSQ and ATLAS [^ experiments have excluded gluinos and squarks up to 1.4 TeV for 
rrig ~ rriq and gluinos up to ~ 0.8 TeV, for the case of rrig 3> m-glQ. 

Many experimental analyses have been performed within the framework of the mSUGRA 
(or CMSSM) model, which assumes a common mass parameter mo (renormalized at the 
GUT scale) for all scalars, and likewise a common mass parameter mi/2 for the gauginos. 
The physical spectrum- obtained by renormalization group (RG) running of soft mass pa- 
rameters from Mgut to M^eak" is characterized by a squark mass spectrum with rriq ~ rrig 
for low mo ~ ^711/2 or m,q 3> m,g for large values of mo- Despite the fact that superpotential 
Yukawa interactions reduce third generation squark masses relative to those of first/second 
generation squarks, third generation squarks nonetheless frequently have masses ~ 1 TeV, 
putting them in conflict with electroweak fine-tuning constraints (discussed below). This 
has led some physicists to question whether weak scale SUSY indeed stabilizes the gauge 
hierarchy, given the constraints from the LHC.^ 

We emphasize that while the various squarks end up being nearly degenerate within 
the mSUGRA model, the limit on m,q quoted above arises mainly from the production 
and decay of first- generation squarks. Only these squarks can be pair produced from the 
valence u and d quark content of the colliding protons. As a result, their production cross 
section falls off much less rapidly with increasing squark mass than the corresponding 
cross section for the production of second and third generation squarks: thus, the ATLAS 
and CMS limits, rUg ~ 1 TeV, apply to first generation squarks, while second and third 
generation squarks may be much lighter without being in conflict with either LHC data or 
with the notion of superpartners as the new physics that stabilizes the weak scale. 

1.2 Impact of LHC Higgs searches 

Recent results from LHC Higgs searches find tantalizing hints for a Standard Model (SM)- 
like Higgs boson of mass rufi ~ 125 GeV (although at present values of ruh ~ 120 GeV 
are also possible). Such a large value of rrih is difficult to realize in models such as mini- 
mal anomaly mediation (mAMSB) or minimal gauge mediation (mGMSB)0, ^j unless all 
sparticle masses are in the 10-20 TeV range, in severe conflict with electroweak fine-tuning 
constraints. Meanwhile, gravity mediation (SUGRA) remains a possible venue for com- 
munication of SUSY breaking since, unlike in mAMSB and mCMSB models, the scalar 
trilinear soft SUSY breaking coupling Ao is an independent parameter, and can be chosen 
to be large as seems to be required by such large values of m/iP]. 



""^By adjusting the trilinear soft breaking parameter Ao to certain values, the ti mass may be dialed to 
sub- TeV values. However, the remaining third generation squark masses typically remain at large values 
and still in possible conflict with fine-tuning constraints. 



SUSY models based on gravity-mediation are, however, expected generically to give 
rise to large FCNC and CP-violating processes [p!o| since there is no mechanism to enforce 
the required generational universality jl|, ^ ^ or alignment of fermion and sfermion mass 
matrices needed to reduce flavour-changing processes to an acceptable level. ■^ Indeed, the 
SUSY flavor and CP problem endemic to gravity-mediation has served as motivation for the 
construction of AMSB and GMSB models, since SUSY sources of FCNC and CP-violation 
are automatically suppressed in these models. 

An alternative solution to the SUSY flavor and CP problems arises by decoupling: 
allowing for first and second generation squark and slepton masses to be in the 10-50 TeV 
range. Third generation squark masses, which directly enter into electroweak fine-tuning, 
or "naturalness" considerations (see below), may be much lighter since flavor and CP 



constraints are relatively mild for third generation particles] 11]. Supersymmetric mod- 
els containing a split spectrum ~ sub- TeV third generation squarks but with multi-TeV 
first/second generation squarks - have been advocated for some time under the label of 
effective SUSY^^, or ESUSY||l^. Indeed, the non-observation of squarks and gluinos in 
the LHC data sample could be a hint pointing in this direction. 

1.3 Naturalness constraints 

It is well known [^ that at tree-level the magnitude of the higgsino mass parameter n 
is determined in terms of (1) the weak scale soft SUSY breaking (SSB) mass parameters 
rn^j^ and m?^ that appear in the Higgs sector scalar potential, (2) the ratio tan/3 = ^, 
and (3) the observed value of the Z-boson mass. Including radiative corrections via the 
effective potential method, this relation gets modified to: 

2^^- (tan2/3-l) ^ ' ^^■^> 

Here, S„ and S^ arise from radiative corrections [O], and are given in the 1-loop approxi- 
mation to the Higgs effective potential by 

1 dAV 



Vu,d ClJ^u,d 

where AV is the one-loop correction to the tree-level potential, and the derivative is eval- 
uated in the physical vaccuum: i.e. the fields are set to their vacuum expectation values 
after evaluating the derivative. 

It is reasonable to say that the theory naturally yields the correct value of Mz if the 



individual terms on the right hand side of Eq. (1.1) are comparable in magnitude so that 
the observed value of Mz is obtained without resorting to large cancellations. Indeed 
this is why |//| has been suggested as a measure of naturalness [|l^] , with theories where 
^2 ^ ]^2 ]-,gj]2g lYiQ "most natural" . This relationship must be accepted with some latitude, 
since values of ff^ ~ (100 GeV)^ are phenomenologically excluded by the LEP2 limit that 



^In mSUGRA, the SUSY GIM mechanism is imposed by simply assuming universaUty of scalar masses 
at the high scale, usually taken to be A/gut- 



rriTT^ > 103.5 GeV. Of course, there is nothing special about fi^ and the same considerations 
apply equally to all the terms, including those involving the radiative corrections. 

In the following, we will somewhat arbitrarily require that each individual term in 
( |1 . 1| ) is bounded by about (200 GeV)^. Similar considerations have recently been adopted 
by several other groups[16, ll^, |l8|, |l9|, pO|.^ In distinction with other works, our focus is 



on the expected sparticle mass spectra and collider and dark matter phenomenology of 
Natural SUSY models with parameters defined at a high scale (taken to be Mgut) which 
lead to weak scale parameters that are natural in the sense that we have just described. 
The largest contributions to T,u^d in Eq. (1.1) arise from superpotential Yukawa interac- 



tions of third generation squarks involving the top quark Yukawa coupling. The dominant 
contribution to these quantities is given by 



m 
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X m? (ln{m^2/Q^) - l) 



and so grows quadratically with the top squark masses. Clearly, the top squark (and by 
SU{2) gauge symmetry, also b^) masses must then be bounded above by the naturalness 
conditions. In Ref. |2C], it has been shown that requiring E^ ~ 2-^1 leads to ?7i^. ~ 
500 GeV. Scaling this to allow /i values up to 150 GeV (200 GeV) leads to a corresponding 
bound rrij:, ~ 1 TeV (1.5 TeV), which of course also applies to the heavier top squark.^ 
In other words, from this perspective, models with fi ~ 150 — 200 GeV and top squarks 
at the TeV scale or below are completely natural. In this connection, it is perhaps worth 
remarking that since 

ml ~ 2//2 + mj^^ + mj^^ + E„ + E^ , (1.2) 

for moderate to large values of tan/3, the heavier Higgs scalars can naturally be at the 
several- TeV scale because of the appearance of tan^ /3 — 1 in the denominator of Eq. (|l.lj). 
Notice, however, that the bound of (200 GeV)^ on each term in Eq. (|1 . 1| ) translates to an 
upper bound 

mA '-^ \mjj^\^ r^ \fi\ tan (3 . (1.3) 

Such a constraint could prove theoretically significant in considerations of high scale models 
with special properties such as models with unified Yukawa couplings at the GUT scale pl|. 



Our discussion up to this point shows that SUSY models with |/i| ~ 150 — 200 GeV and 
top squark masses (and if tan /3 is very large, also bottom squark masses) below 1-1.5 TeV 
are perfectly natural. There will also be corresponding constraints on other sparticles such 
as electro- weak charginos and neutralinos that also directly couple to the Higgs sector, 
but since these couplings are smaller than ft and because there are no colour factors, the 
constraints will be correspondingly weaker. Sparticles such as first and second generation 
squarks and sleptons that have no direct/significant couplings to the Higgs sector are 



^These analyses differ in detail on the restrictions on each term, and even whether a common constraint 
is applied to each term. For this reason, some of the constraints that have been obtained by these analyses 
are stronger than the ones we obtain in this paper. 

*For very large values of tan/? where the bottom quark Yukawa becomes comparable to ft, similar 
considerations would apply to bottom squarks. 



constrained only via two-loop effects and can easily be in the 10 - 50 TeV range. An 
important exception would be the gluino, since radiative corrections to the top squark 
mass are proportional to ?tt.^|17|. Using (5m? ~ a^'^i ^ ^og and setting logs to be order 

unity, we see that nig ~ Srriq. For top squarks to remain below the 1.5 TeV range, the gluino 
must be lighter than about 4 TeV. In models with electroweak gaugino mass unification, 
electroweak-inos would then automatically not destroy naturalness. 

1.4 Natural SUSY 

These considerations suggest that the region of SUSY parameter space where 

• 1^1 ~ 150 - 200 GeV, 

• third generation squarks m^ , m^ ~ 1 — 1.5 TeV, 

• rrig ~ 3 — 4 TeV and SSB electroweak-ino masses smaller than 1-2 TeV 

• niA ~ l/^l tan/3, 

• rriq^^, m^ ~ 10 — 50 TeV, 

may, from naturalness and flavor /CP considerations, merit a dedicated study. The first 
and second generation squarks and sleptons - lying in the 10-50 TeV range - provide 
a decoupling solution to the SUSY fiavour problem, the SUSY CP problem and to the 
problem of too-rapid-proton decay. We remark here that if SUSY breaking arises from 
supergravity breaking in a hidden sector, then the gravitino mass m.3/2 sets the scale for 
the largest of the SSB terms, and we would also expect m^/2 ~ 10 — 50 TeV: such a high 
value of m.3/2 also provides a solution to the gravitino problem 1 22, ^]. The heavier Higgs 



bosons may easily be in the several- TeV range for moderate to large values of tan/3. 

SUSY models with the above generic spectra have been dubbed "Natural SUSY"|18| 
and are a more restrictive case of effective SUSY models because we further restrict \fi\ ~ 
150 — 200 GeV. This usually gives rise to a higgsino-like lightest neutralino Zi. 

1.5 Dark matter in natural SUSY 

In fact, a problem with effective SUSY models with a bino-like Zi arises in that a vast 
overabundance of neutralino cold dark matter (CDM) is expected [p^], typically 2-4 orders 
of magnitude above the WMAP-measured value of Q.cDKih'^ ~ 0.11 unless weak scale 



parameters happen to be in special parameter space regions. It has been suggested in |13| 
that if the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism - which 
introduces a supermultiplet containing spin-zero axion and saxion fields, along with a 
spin-i axino - then neutralinos might decay to a light axino LSP via Zi — )• 07. Since 
each neutralino converts to one axino, the decay-produced axino abundance is given by 
^d'^^ = rria/m^ ^^ h? ■ For rria in the MeV range, the suppression factor is ~ 10^^ — 10"^, 
bringing the DM abundance into accord with measurement. 

However, typically in gravity mediation models the axino mass is expected to be around 
the TeV scale[^, ^, with Zi remaining as LSP. In fact, in the PQ-augmented SUSY 



model, one then expects thermal axino production (TP) in the early universe, followed by 
late-time a — )■ Z\^ decays, so the dark matter overabundance is made even worse. In the 
case of natural SUSY, the higgsino-like Z\ leads typically to a thermal underabundance of 
neutralino CDM. But now TP axinos followed by their decay to neutralinos can augment 
this abundance 1 2(:, ^, 28 1, while any remaining underabundance can be filled by axions 



produced via vacuum mis-alignment (coherent oscillations). Thus, in this case we might 
expect the CDM to consist of a higgsino-axion admixture. Which of these two particles 
dominates the DM abundance depends on specific choices of PQ parameters and on the 
value of re-heating temperature Tpt after inflation. 

2. Parameter space and mass spectra for Natural SUSY 

Since the introduction of softly broken SUSY into the Standard Model (leading to the 
Minimal Supersymmetric Standard Model, or MSSM) leads to a theory with stable mass 
hierarchies, it is natural to assume the MSSM is the low energy effective theory arising 
from an underlying SUSY Grand Unified Theory (GUT) which is broken at some high 
energy scale, taken here for definiteness to be Mqut — 2 x 10^^ GeV. Indeed, the MSSM 
(or MSSM plus gauge singlets and/or additional complete SU{5) multiplets) receives some 
indirect support from experiment in that 1. the measured weak scale gauge couplings unify 
nearly to a point at Mqut under MSSM renormalization group (RG) evolution and 2. the 
MSSM electroweak symmetry is broken radiatively due to the large top quark Yukawa 
coupling, consistent with the measured value of mt- 

Motivated by these SUSY success stories, the interesting question arises as to whether 
the natural SUSY sparticle mass spectrum can be consistently generated from a model 
with parameters defined at the high scale Q = Mqut- To implement a low value of |/i| ~ 
150 — 200 GeV, we will adopt the 2-parameter non-universal Higgs model (NUHM2)[pg|], 
wherein weak scale values of /i and rriA may be used as inputs in lieu of GUT scale values 
of m\j and rn^j . To generate the split first/second versus third generation scalar mass 
hierarchy, we will adopt a common GUT scale mass mo (3) for third generation scalars, 
and a common GUT scale mass mo (1,2) for the first/second generation scalars. The 
intra-generational mass universality is well-motivated by 5*0(10) GUT symmetry, since all 
matter multiplets of a single generation live in the 16-dimensional spinor rep of 5*0(10). 
We can also allow some degree of non-universality between m-o(l) and m,o(2) so long as 
both lie in the tens of TeV regime, and provide a decoupling solutions to SUSY FCNC 
and OP- violating processes (for constraints from FCNC processes, see Ref. [^0|). For 
convenience, we will take them as degenerate. 

To allow for a light third generation, we adopt different GUT scale values for the scalar 
mass parameter of the first two generations and the third generation. In the spirit of SUSY 
GUT theories, we will assume gaugino mass unification to a common gaugino mass mi 12, 
and assume a universal trilinear scalar coupling ^o a-t the GUT scale. The sparticle mass 
spectrum together with sparticle couplings is then determined by the parameter set,^ 

mo(l,2), mo(3), mi/2, ^o, tan/3, /i, tua ■ (2.1) 



^We could have allowed a D-term contribution that would arise when the additional (/(I) that is in 



We take mt = 173.3 GeV. 

Our goal in this Section is to search for weak scale spectra that are natural in the sense 
defined above within this framework, to delineate regions of parameter space consistent with 
low energy constraints, and to study their implications for SUSY searches at the LHC. 

For our calculations, we adopt the Isajet 7.82 [Q SUSY spectrum generator Isasugra|32] 



Isasugra begins the calculation of the sparticle mass spectrum with input DR gauge cou- 
plings and /b, /r Yukawa couplings at the scale Q = Mz {ft running begins at Q = mt) 
and evolves the 6 couplings up in energy to scale Q = Mqut (defined as the value Q 
where gi = (72) using two-loop RGEs. We do not strictly enforce the unification condition 
93 = 91 = 92 at Mgut, since a few percent deviation from unification can be attributed to 



unknown GUT-scale threshold corrections [33|. At Q = Mgut, we introduce the SSB pa- 



rameters in ( |2.1| ) as boundary conditions, and evolve the set of 26 coupled MSSM RGEs |34| 
back down in scale to Q = Mz- Full two-loop MSSM RGEs are used for soft term evo- 
lution, while the gauge and Yukawa coupling evolution includes threshold effects in the 
one- loop beta- functions, so the gauge and Yukawa couplings transition smoothly from the 
MSSM to SM effective theories as different mass thresholds are passed. In Isasugra, the 
values of SSB terms which mix are frozen out at the scale Q = MgusY = J'^^t ^-t ; while 
non-mixing SSB terms are frozen out at their own mass scale |£^]- The scalar potential 
is minimized using the RG-improved one-loop MSSM effective potential evaluated at an 
optimized scale Q = MsusY which accounts for leading two-loop effects [^]. Once the 
tree- level sparticle mass spectrum is computed, full one-loop radiative corrections are cal- 
culated for all sparticle and Higgs boson masses, including complete one-loop weak scale 
threshold corrections for the top, bottom and tau masses at scale Q = MsusY ^^- Since 
the GUT scale Yukawa couplings are modified by the threshold corrections, the Isajet RGE 
solution must be imposed iteratively with successive up-down running until a convergent 
sparticle mass solution is found. Since Isasugra uses a "tower of effective theories" approach 
to RG evolution, we expect a more accurate evaluation of the sparticle mass spectrum for 
models with split spectra (this procedure sums the logarithms of potentially large ratios of 
sparticle masses) than with programs which make an all-at-once transition from the MSSM 
to SM effective theories. 

We searched for Natural SUSY solutions in the above parameter space by first fixing 
fi = 150 GeV, and then performing a (linearly weighted) random scan over the remaining 
parameters in the following ranges: 

mo(l,2) : 5-50 TeV, (2.2) 

mo(3) : 0-5 TeV, (2.3) 

mi/2 : - 5 TeV, (2.4) 

-4< Ao/moiS) <4, (2.5) 



50(10) but not in the SU{5) subgroup is spontaneously broken. However, because we allow jj, and tua as 
inputs, this would, however, have no impact on the allowed range of nih, and so would only affect sparticle 
masses. In order to avoid a time-consuming scan over yet one more parameter which would probably not 
qualitatively alter the main results presented below, we have not included this term in our parametrization 
of the model. 



ruA ■■ 0.15 - 2 TeV, (2.6) 

tan/3 : 1-60. (2.7) 

We require of our solutions that (1) electroweak symmetry be radiatively broken (REWSB), 
(2) the neutrahno Zi is the hghtest MSSM particle, (3) the light chargino mass obeys the 



rather model independent LEP2 limit that m^ > 103.5 GeV|37] and (4) that rrig < 4 
TeV, in accord with our naturalness criterion detailed above. 

The results of our scan are plotted in Fig. ||. On the y-axis, we plot the average third 
generation squark mass 

rng{3) = {mi^ + m^^ + mjJ/3 (2.8) 

while the x-axis lists the particular parameter. Blue points have mg(3) < 1.5 TeV, green 
points have rrTg(3) < 1 TeV and red points have rnq{3) < 0.5 TeV. 

In frame a), we see that we can generate solutions with fnq{3) lower than 0.5 TeV, 
but only for values of mo{l,2) ~ 18 TeV. For heavier values of mo(l,2), it is well known 
that two-loop RGE effects tend to push third generation squark masses into the tachyonic 



range pq, 39, |T^, which here would correspond to color breaking minima in the scalar 
potential. On the other hand, requiring mg(3) < 1 (1.5) TeV allows for ?tt,o(1; 2) as high as 
~ 25 TeV- enough to suppress FCNCs except in the case of very large flavor violating soft 



terms [38 1 . In frame 6), we plot the required value of mo (3) to give rise to sub-TeV average 



squark masses: here, values of mo (3) < 2 (5) TeV are required to generate solutions with 
rnq{3) < 0.5 (1) TeV. Frame c) shows the value of mi/2 required for natural SUSY models. 
A value of mi/2 < 1-4 TeV is required for mg{3) < 0.5 TeV, while mi,2 ~ 1.7 TeV because 
we impose nig ~ 4 TeV. In frame d), we see that rnq(3) < 0.5 TeV can only be achieved 
for ^0 ~ 0, while rnq(3) < 1 TeV is allowed for Aq > — 2mo(3), i.e. Aq cannot be large, 
negative. In frame e), we find that rnq{3) < 0.5 is allowed for tan/3 < 50, while rnq{3) < 1 
TeV can be achieved for any tan /3 from ~ 2 - 60. Finally, frame f) shows that solutions 
with rnq{3) < 0.5 TeV can be found for any value of niA '■ 0.15 — 2 TeV. 

In Fig. |2|, we show the value of ruh which is generated in NS models versus various 
SUSY parameters. In frame a), we see that the red points with mq{3) < 0.5 TeV populate 
the range rrih ~ 105 — 120 GeV, while mh values, as obtained using Isajet, up to 123 GeV 
(124 GeV) can be readily accommodated for mq;(3) up to 1 TeV (1.5 TeV). This should 
be compared with 115.5-131 GeV (114-127 GeV), the range of light Higgs boson masses 
currently allowed by the ATLAS (CMS) data [||, H at the 95%CL. These experiments 
also report a small excess of a signal at ruh ~ 125 GeV. For the smallest range of rnq{3) in 
the figure, it might appear that one would be hard pressed to accommodate the LHC hint of 
a 124-126 GeV light Higgs scalar. Of course, here one must keep in mind that Atlas/CMS 
may really be seeing a Higgs scalar with mass closer to 124 GeV or that mt may be slightly 
larger than 173.2 GeV as assumed in our calculation of the radiative correction. Combining 
this with a ~ 3 GeV error anticipated in the Isasugra calculation of m/j and it becomes 
apparent that values of nih ~ 120 — 121 GeV may be consistent with the Atlas/CMS /i(125) 
hint even for small values of mq;(3). The largest values of rrih are obtained for tan/3 ~ 10. 
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Figure 1: The value of mq{3) versus various SUSY parameters with /i = 150 GeV. The dots are 
colour-coded by the range of TOq(3):< 0.5 TeV (red); 0.5-1 TeV (green); 1-1.5 TeV (blue). 



We have already seen in Fig. ^d) - and also here in Fig. ^d)- that only Aq > 
values lead to rnq{3) < 0.5 TeV, while in Ref. [^, ^ it is found that the largest values 
of nifi are found for Aq ~ —2mo. As we allow increasing values of mg(3) consistent with 
our naturalness conditions, we see that values of ^o ~ — 2?n,o(3) become allowed, and 
consequently higher values of mh can be accommodated. This is the case of maximal 
mixing in the top squark sector, which leads to maximal nih values p2|. 



The value of BF{b — t- S7) should be rather tightly constraining for models of natu- 
ral SUSY, since there may be several light third generation squarks, and not too heavy 
charginos, and since the main non-standard contributions to the decay rate come from 

to 



top-squark-chargino loops [ 43 1. Here, we implement the Isatools subroutine IsaBSG| 
compute the branching fraction, which is listed in Figure ^ versus mq{3). These values are 




H 



120- 



I*" 

115- ih, 

- * • 

110 '.. . 

105 -.\r ■ 






' ^ ' ■vS**'^^^^''^'^''-^-'^?'^**'*^^^^^^^^ ■ ^ - '^ ■• ■• ■"• ■ ''■' ■ ■ 



■••.',:\'J, ■■,'':.:'. .. .• • 

: ■..■-■;#; ■v•:^■. • ■ • 



~130 



E 



m„(3) (TeV) 



'"■''■■'^^^^^ 



■ B- *> .V,.. 

•.,< "^^ ,J,I 






3 4 

Ao/m„(3) 





: 


125 


- 














- 






120'^ 








' . 


-i 








^ f 


115 >• 


^ 






•V 


110 ; 




, » 


/ 




105 


' 1 • 


r 


1 , 

. r> . *. .1 .-. . 1 


1 . 11 1 . 11 '1 


< < 1 < < < 



tanp 



0.2 0.4 0.6 0.8 



1.2 1.4 1.6 1.8 2 

m, (TeV) 



Figure 2: The value of mh versus various GUT scale SUSY parameters. Here, and in subsequent 
figures, the colour-coding is as in Fig. ^ 



to be compared with the measured value of BF{b — )■ 57) = (3.55 ± 0.26) x 10~^ from Ref. 
IQ. Indeed, as can be seen, large SUSY loop contributions cause the branching fraction to 
vary over a wide range: (0 — 9) x 10~^, so that many solutions would be rejected. Nonethe- 
less, many other solutions do remain within the itSo" band (which is shown), where the 
various loop contributions may cancel one against another to yield consistency with the 
measured value. 

In addition, the well-known (g — 2)^ anomaly has been reported as a roughly 3cr 



deviation from the SM value: Aa^ = (28.7 ± 8.0) x 10 



-10 r 



In Natural SUSY, since the 



/ii^2 and Dfj^ masses are in the multi-TeV range, only a tiny non-standard contribution to 
the {g — 2)^ anomaly is expected, and alternative explanations for this anomaly would have 
to be sought. 
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Figure 3: Predicted values of the branching fraction for b —> sj vs. mq{3). We also show the 
experimentally determined central value ±3cr band for the BF{b — > 57). 



3. Benchmark points, slopes and planes 

In this section, we list some representative natural SUSY benchmark points, slopes and 
planes which could be used for LHC analyses. In Table [l|, we show three such points, NSl 
with mg{3) = 666 GeV, NS2 with mg(3) = 595 GeV and NS3 with mg{3) = 1343.7 GeV. 
For all points, we fix /i = 150 GeV, with large w-1/2 so that the higgsino-like chargino 
and the two lightest higgsino-like neutralinos have masses ~ 150 GeV and are the lightest 
sparticles. The light Higgs masses ruh ~ 121 GeV for the first two points, and so are 
low but as discussed above not incompatible with the recent hint for ruh ~ 125 GeV. The 
third point NS3 allows rrih = 123.5 GeV but at the expense of rather large rnq{3), and 
mr marginally beyond our naturalness requirement. For all these points the gluino mass 
is around 3 TeV and first and second generation squarks are completely beyond the reach 
of the LHC. 

For point NSl, the light top squark ti is next-lightest SUSY particle after the three 
higgsino-like states; it has mass m^ = 301.4 GeV and may be accessible to LHC top 
squark searches. The 61 and t2 come in at 788 and 909 GeV, respectively. Both staus 
and the tau sneutrino are relatively light and might be accessible at a future TeV-scale 
lepton-anti-lepton collider. 

Point NS2 has a light bottom squark with m^ = 497.3 GeV as next-lightest after the 
higgsinos. The ti is slightly heavier at 572 GeV. This point has heavier tau sleptons which 
would not be accessible to any planned lepton colliders. Point NS3 with rather heavy third 
generation squarks and sleptons would be very challenging to see at LHC although the 
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parameter 


NSl 


NS2 


NS3 


mo(l,2) 




13363.3 


19542.2 


7094.3 


mo(3) 




761.1 


2430.6 


890.7 


"1-1/2 




1380.2 


1549.3 


1202.6 


^0 




-167.0 


873.2 


-2196.2 


tan/3 




22.9 


22.1 


19.4 


f^ 




150 


150 


150 


rriA 




1545.6 


1652.7 


410.1 


Vflg 




3272.2 


3696.8 


2809.3 


rriuL 




13591.1 


19736.2 


7432.9 


mun 




13599.3 


19762.6 


7433.4 


ruej. 




13366.1 


19537.2 


7086.9 


i^h 




301.4 


572.0 


812.5 


^h 




909.2 


715.4 


1623.2 


^h 




788.1 


497.3 


1595.5 


^b2 




1256.2 


1723.8 


1966.7 


TTtfi 




430.9 


2084.7 


652.2 


mf^ 




532.9 


2189.1 


1065.5 


rriy^ 




402.3 


2061.8 


1052.1 


^w^ 




1180.2 


1341.2 


1013.9 


"^w. 




155.9 


156.1 


156.2 


^z^ 




1181.3 


1340.4 


1020.0 


'"fa 




615.3 


698.8 


532.6 


"^^2 




156.3 


156.2 


157.0 


^Z^ 




148.4 


149.2 


147.4 


rrih 




121.3 


121.1 


123.5 






0.007 


0.006 


0.007 


BF{b^ 


Si) 


2.8 X 10-4 


3.6 X 10-4 


2.8 X 10-4 


a^\Z^p) (pb) 


5.5 X 10-9 


1.8 X 10-9 


9.8 X 10-9 


a^^iZ^p) (pb) 


3.9 X 10-5 


2.9 X 10-5 


5.7 X 10-5 


(o-w)|i;^0 


(cm^/sec) 


3.0 X 10-25 


3.1 X 10-25 


3.0 X 10-25 



Table 1: Input parameters and masses in GeV units for three Natural SUSY benchmark points, 
with n — 150 GeV. Also shown are the values of several non-accelerator observables. 



spectrum of light higgsinos should be accessible to a linear e^e collider. 

In Fig. El we convert benchmark point NSl into a benchmark slope by retaining 
all parameters as in Table IT], except allowing mo (3) to vary. For mo (3) much below 700 
GeV, we generate spectra with tachyonic stops. Some gaps occur in the plot where no 
convergent RGE solution is found. These gaps can be filled in by increasing the number 
of iterations in Isasugra RGE running beyond the default value of 25. In the Figure, we 
plot all four third generation squark masses versus mo (3), which gives a rising spectrum for 
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most third generation squarks except the light top squark which reaches a minimal mass 
at 7TT-o(3) — 840 GeV, where rrif < m^ so that the ti is lightest MSSM particle. This 
point gives maximal mixing in the top squark sector, and a minimal value for mr . 
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Figure 4: Plot of third generation squark masses together with the lightest neutralino mass versus 
variation in 7tt,o(3), with other parameters fixed as for the benchmark point NSf . 

In Fig. |5|, we convert benchmark point NSl into a benchmark plane, where we plot 
contours of light top-squark mass mi as a function of mo(3) vs. fi variation. The unshaded 
region gives rise to tachyonic squarks. Points with valid solutions are labeled as black 
dots; the gaps again require iterations beyond the default value of 25. The color coding 
extrapolates the generated value of m^ , which again reaches a minimum of below 200 GeV 
at mo (3) ~ 830 GeV. 

4. LHC signals for Natural SUSY 



We begin by noting that since /x ~ 200 GeV, we expect a spectrum of light, higgsino-like 
Wi, Z\ and Z^ with mass ~ /i and small mass gaps m^ — m?^ ^ m^ — m;^ ~ 10 — 20 



m^ ~ m^ 



GeV. Models with low fi parameter have been considered previously in Refs. [46| and |47|. 
In Ref. 1 47], production cross sections for chargino pair production, chargino- neutralino 
production and neutralino pair production were presented. The -ino pair production cross 
sections tend to be in the 50-500 fb range. The decays Wi — )■ Ziff and Z2 — )■ Ziff 
(where / collectively stands for light SM fermions) are dominated by W* and Z* exchange 
respectively. However, since the mass gaps Wi — Z\ and Z2 — Z\ are so small, there is 
only a small visible energy release in the decays, making the visible portion of the final 
state very soft and difficult to extract above SM backgrounds. In fact, models with low 
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Figure 5: The top squark mass m^ in tlie fi vs. too(3) plane, witli otlier parameters as for 
benchmaric point NSl. 

fj, and concomitantly light higgsinos but other sparticles at the multi-TeV scale have been 



dubbed "hidden SUSY" in Ref. [47| because distinctive SUSY signals at LHC are extremely 
difficult to extract above background. 

A key feature of Natural SUSY models is that they necessarily feature three and 
possibly four relatively light third generation squarks. While simplified models tend to 
focus on the signal from a single production mechanism, often assuming one dominant 
decay channel, generally speaking in natural SUSY we expect several third generation 
squarks to contribute to new physics signal rates. Moreover, these squarks will typically 
have decays to all three light higgsino-like states, and possibly also to other decay channels. 
While the lightest of these squarks will have the largest production cross sections, because 
of the larger mass gaps together with their more complex cascade decays, production of 
heavier third generation squarks may also yield observable signals. 

In Fig. P, we list the pp — )• titiX production cross section calculated in NLO QCD 
using Prospino|4^. We show results for LHC with ^/s = 7, 8 and 14 TeV center-of-mass 
energy.^ The bibi (for i = 1,2) and ^2^2 cross sections are essentially identical to those 
shown by making an appropriate mass substitution, since almost all the production cross 
section comes from light quark qq and gg fusion in the initial state. 

In Table ^, we show the various third generation squark pair production cross sections 
at LHC8 and branching fractions for benchmark points NSl, NS2 and NS3 from Table 



®More precisely, because Prospino only allows a selection of Tevatron, LHC7 or LHC14 ~ but not LHC8 
- we have obtained the cross section for LHC8 by scaling the corresponding Isajet cross section by the ratio 
of Prospino to Isajet cross sections for LHC7. 
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Figure 6: Top squark pair production cross sections at LHC7, LHC8 and LHC14 versus m^ . Witli 
the appropriate mass substitution these curves also represent the cross sections for pair production 
of ^2, bi or 62 squarks. 



|l[ Point NSl is by far dominated by iiii production at LHC8 with a cross section of 
~ 2 pb. The ti then decays to bWi at ~ 100% branching fraction. This model would 
be well-described by a simplified model analysis, where the final state Wi is essentially 
regarded as missing- £^r due to its soft decay products. Thus, the signature would be a 
pair of acollinear 5-jets together with E^^^^ and no other transverse activity except from 
QCD radiation. The dominant SM physics background would be from Zbb production, 
with Z — )• uv'. At LHC8, there is also a 4 fb cross section from bibi production followed 
by bi — )• Wti, giving rise to bbW~^W~ + £'™^^ events, albeit at low rates. These would be 
subject to a daunting background from tt production. We mention that at LHC14, t2-pair 
production which has a cross section of ~ 20 fb could lead to a handful of spectacular 
^2^2 — ^ ZZtiti — )• ZZbb + EJ^^^^ events where the Zs might be identified via their leptonic 
decays. 

For the point NS2, pp — t- bibi production is dominant at a ~ 80 fb, although titi is 
also produced at ~ 30 fb. In this case, the bi decays dominantly via bi — t- tWi giving 
rise to a it -|- E^^^^ signature at LHC. The decays bi — >• bZi and 61 — ?■ 6Z2 also occur at 
~ 10% level. The ii decay modes are spread somewhat evenly between bWi, tZ\ and tZ2 
final states. Again, the small mass gap between the W\l Z2 and the LSP implies that the 
chargino and the neutralino daughters are essentially invisible. By combining all modes, 
the most lucrative signature channels consist of 66 -|- E^^^^ and ti + E^^^^ events. The 
heavier t2 decay modes are spread among many more possibilities, including decays to W 
and h bosons in the final state; a handful of novel events may be obtained at LHC8, but 
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parameter 


NSl 


NS2 


NS3 


a{iiii) 


2000 fb 


30 fb 


2fb 


BF{ii -^ bWi) 


1.0 


0.25 


0.62 


BF{ii -^ tZi) 


- 


0.42 


0.08 


BF{ii -^ tZ2) 


- 


0.33 


0.30 


CT{hlbl) 


4fb 


80 fb 


0.00013 fb 


BF{bi -^ bZi) 


0.01 


0.10 


0.01 


BF(bi -^ bZ2) 


0.01 


0.09 


0.01 


BF(bi -^ tWi) 


0.09 


0.81 


0.04 


BF{bi -^ Wti) 


0.89 


- 


0.94 


(Jihh) 


Ifb 


6fb 


0.00011 fb 


BF{i2 -^ bWi) 


0.09 


0.29 


0.05 


BF{i2 -^ Zti) 


0.70 


0.01 


0.39 


BF{i2 -^ hii) 


0.01 


0.23 


0.25 


BF{i2 -^ Wbi) 


0.03 


0.16 


0.26 


BF{t2 -^ tZi) 


0.09 


0.13 


0.03 


BF{i2 -^ tZ2) 


0.08 


0.16 


0.02 


crihh) 


0.05 fb 


0.0001 fb 


0.00004 fb 


BF(b2 -^ bZi) 


0.22 


0.23 


0.01 


BF{b2 ^ bZ2) 


0.22 


0.22 


0.01 


BF{b2 -^ bZs) 


0.07 


0.08 


- 


BF(b2 -^ tWi) 


0.42 


0.44 


0.02 


BF(b2 -^ Wii) 


0.03 


0.01 


- 


BF(b2 -^ hbi) 


0.03 


0.02 


- 


BF{b2 -^ Hbi) 


- 


- 


0.23 


BF{b2 -^ All) 


- 


- 


0.23 


BF(b2 -^ H-ii) 


- 


- 


0.41 


BF(b2 -^ H-t2) 


- 


- 


0.08 



Table 2: Production cross sections at LHC8 and branching fractions for third generation squark 
production for the Natural SUSY benchmark points from Table 1. 



more likely at LHC14. The 62 state appears likely undetectable even at LHC14. 

For the benchmark point NS3, the detection of third generation squarks at LHC8 
appears to be very difficult on account of the very low cross sections. Even at LHC14, the 
cross section for titi production is just 50 fb, and the fact that the chargino and neutralino 
daughters are (nearly) invisible will make identification of the acollinear ti and bb events 
from this quite challenging. Production of ^2^2 and bibi at LHC14 occurs at a fraction of a 
fb level, though the interesting topologies that include Z and h production from ^2 cacsade 
decays may be accessible at super-LHC luminosities. 

We would also like to assess the prospects for discovering the gluino of the natural 
SUSY framework at the LHC. A plot of rrig vs. mq;(3) is shown in Fig. ^ for scan points 
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fulfilling the BF{b — t- 57) constraint and also rrih > 115 GeV. From the plot, we see that 
while models with rrig ~ 1 TeV can be readily obtained for rnq{3) ~ 1 — 1.5 TeV, a tighter 
restriction of rnq{3) < 0.5 (1 TeV) typically limits rrig ~ 2 TeV (1 TeV). In the case of 
models with multi-TeV squarks, the LHC8 reach (which should be close to LHC7 reach pS|]) 
for ~ 20 fb^^ fb of integrated luminosity extends to rrig ~ 1 TeV. The LHC14 reach [^] for 
100 fb~^ extends to ~ rUg ~ 1.8 TeV. These studies have been done within the mSUGRA 
model, and for LHC14 without 6-jet tagging which should enhance the SUSY signal in 
Natural SUSY models. The increased reach in the gluino mass is projected to be up to 



20%, depending on the details of the particle spectrum |5l|]. We conclude that while some 
models with large fnq{3) > 1 TeV may be accessible to LHC gluino searches,^ there remain 
many models (expecially for low rnq{3) < 1 TeV) where gluino pair production will be 
beyond even the LHC14 reach. 
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Figure 7: Value of rrig vs. fnq{3) from Natural SUSY models which obey BF{b — > 57) at 3cr and 
rrih > 115 GeV. 



Lastly, motivated by the bound on rriA presented in Sec. |L^, we plot rriA vs mq{3) in 
Fig. ^ for natural SUSY points with rrig < 4 TeV, rrih > 115 GeV and which satisfy the 
B{b — )• 57) constraint. The color coding is as in Fig. 0. We see that niA ~ 500 GeV for 
low rnq{3) values, while rriA can be as low as a few hundred GeV for very large friq{3). 



^In this context, we note that the ATLAS LHC7 hniits |52] from gluino-mediated stop-pair searches do 
not directly apply because these rely on the same-sign dilepton signal where the lepton may arise from 
either the top quark daughter of the gluino or from the chargino daughter of the top squark. In our case, 
we expect leptons from the latter source to be very soft. 
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Figure 8: Value of niA vs. r?T,g(3) in Natural SUSY models with rrig < 4 TeV, mh > 115 GeV and 
which satisfy BF{b — > 57) at 3a. The color coding is as in Fig. |^. 



5. Natural SUSY at a linear e~^e collider 



While Natural SUSY may possibly be difficult to discover at the LHC, it leads to a potential 
bonanza of signals for a linear e~^e~ collider (LC) operating in the 0.3 — 1 TeV range. A LC 
would potentially be a higgsino factory because, as we have emphasized, the ^ parameter 
is necessarily small in this scenario. Indeed, pair production of higgsino-like chargino 
and also neutralino states W^ , Zi and Z2 with sizeable cross sections is inevitable at a 
0.25 — 0.5 TeV machine. Example cross section plots have been shown in Ref. | |47| ] and so 
will not be repeated here. Thus, a natural target for a LC would be the pair production 
reactions e^e~ — )• W^W^ , Z1Z2 and Z2Z2. While the small Wi — Z\ and Z2 — Z\ mass 
gaps are a formidable challenge at LHC (and may also be problematic at a LC), it has 



been shown[53| that with specialized cuts, it should be possible to extract a signal above 
SM background at a LC. The visible energy from these reactions would be low just above 
threshold, but as -\/i increases, the decay products from W\ and Z2 would be boosted to 
higher energies. In addition, the beam polarization would be a strong tool not only for 
distinguishing the signal from VF"*"!^" backgrounds, but also for distinguishing between 
wino-like versus higgsino-like charginos|47]. 

In the case of the NSl benchmark, after the light higgsinos are well studied and the 
CM energy y^ is increased, the next target threshold would be i\i\ production at y^ ~ 
2?7Zj ~ 610 GeV. This would be followed by the tau-snuetrino pair production threshold at 
y^ ~ 810 GeV, with v-,- — )• W^t^ decay. At a little higher energy, t^t^ pair production 
would turn on, followed mainly by f\ — t- Z\t and Z2T decay. For the heavier spectra shown 



17 



in NS2 and NS3, the light higgsino pair production reactions would still be available, 
but CM energies of over 1 TeV would be required to pick up any squark pair production 
reactions. As emphasized above, the accessibility of higgsino-like states is a generic feature 
of Natural SUSY models. 

6. Natural SUSY and direct/indirect WIMP searches 

As noted in Sec. [l|, the higgsino-like neutralinos with masses ~ 100 — 200 GeV expected 
in NS models annihilate very rapidly in the early universe and so yield a thermal relic 
underabundance of CDM. However, the neutralino relic abundance can be boosted to 
match its observed value in models where 

• a PQ solution to the strong CP problem is invoked, and thermally-produced but late- 



decaying axino (and/or saxion) decays augment the SUSY particle production [B^, 28|, 
or 

• there exist late-decaying TeV scale moduli fields with large branching fractions to 
SUSY particles that subsequently decay to the neutralino [p^]. 

We stress that it in the first case it is not necessary that neutralinos saturate the 
observed relic density. Indeed it is easy to select PQ parameters where the converse is 
true: O^ h? stays low while the bulk of dark matter is comprised of axions, or even where 
both axion and neutralino abundances are comparable. In this case, the direct and indirect 
neutralino reach estimates presented below (these have been obtained assuming that neu- 
tralinos saturate the CDM density) would have to be increased by a factor of O.II23/O5 h? . 
It is difficult, but not impossible p8| , to lower the neutralino abundance below its standard 
thermally produced prediction. Thus, we expect the above reach scale factor to be typically 
between 1 and 16, since Q}-~'^^'^''^^h? ~ 0.007 for a higgsino-like relic neutralino with a mass 
~ 150 GeV. 

In Fig. P, we show the spin-independent Zip scattering cross section in pb as obtained 



from IsaReS|55|. Here, and in the remainder of this section we assume that Zi saturates 
the DM density. We plot points versus mg(3), since m^ is fixed typically ~ 150 GeV due 

to our choice of /u = 150 GeV. We actually find that the bulk of points inhabit the ~ 10^® 
pb range. Comparing to the bound from Xe-100[^], we find that a large fraction of these 
points are excluded if the higgsino-like WIMP is essentially all the dark matter. Moreover, 
with this same assumption, a large fraction of surviving points lie within the projected 
reach of Xe-100/2012 run, and certainly within the reach of Xe-lton. 

In Fig. ^, we plot the spin-dependent a {Zip) scattering cross section, this time in 



comparison to current and future IceCube reach|5^, and future COUPP reach[58|. While 
the current IceCube reach excludes a significant portion of points (under the asumption 
of neutralino dominance), the future IceCube and especially COUPP reaches will access 
most of the remaining parameter space. 



Fig. 11 shows the thermally averaged neutralino annihilation cross section times rel- 



ative velocity, evaluated as f — t- 0. This quantity enters linearly into indirect searches for 
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Figure 9: Spin independent pZi scattering cross section versus too(3) for NS models with /i = 150 
GeV. Also shown are the current 90% CL bounds together with projections from the XENONIOO 
(2012 sensitivity) and IceCube (180 day sensitivity) experiments for a 150 GeV WIMP (assuming 
higgsino-like WIMPs saturate the measured dark matter abundance). 



neutralino annihilation in the cosmos into 7s or e"*", p or D searches. For our case, the bulk 
of points inhabit the range (crf)|^,_^o ~ 10^^^ cni'^/sec. The horizontal solid line shows the 
upper limit on the annihilation cross section times velocity for very non-relativisitic dark 
matter in dwarf spheroidal satellite galaxies of the Milky Way annihilating to W boson 



pairs obtained by the Fermi collaboration [59], assuming a ~ 150 GeV WIMP. Models with 
a larger annihilation cross section would have led to a flux of gamma rays not detected 
by the experiment, assuming a Navarro-Frenk- White profile for each dwarf galaxy in the 
analysis. We see that the Fermi bound excludes the bulk of points for our choice of DM 
mass, again assuming higgsinos saturate the DM density. Moreover, this bound changes 
rather slowly with the DM mass, being just a factor of 2 weaker for a WIMP mass of 
300 GeV. Further searches and improvements by the Fermi-LAT Collaboration and/or the 
impending AMS results should provide an incisive probe of the NS framework. 



7. Summary and conclusions 



The Natural SUSY model is defined by distinctive spectra characterized by a low |/x| 
TUfi ~ 150—200 GeV, with a rather light spectrum of third generation squarks m^ , 



t2 



mr ~ 

01 



0.5 — 1.5 TeV to stabilize the electroweak scale. In addition, rrig ~ 4 TeV so that loop cor- 
rections to third generation squark masses are smaller than the squark mass. First/second 
generation sfermions, on the other hand, could be at the tens of TeV scale, thus suppress- 
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Figure 10: Spin dependent pZi scattering cross section versus mo(3) for NS models witlr /i = 150 
GeV. Also shown are current limits and future projections from IceCube and COUPP experiments 
(assuming higgsino-like WIMPs saturate the measured dark matter abundance) . 

ing unwanted flavor-violating and CP-violating processes. Motivated by gauge coupling 
unification, we expect the MSSM, or MSSM plus gauge singlets, to be the effective field 
theory between M„eak and Mqut- In this case, the Natural SUSY mass spectra should 
arise from underlying fundamental parameters that have their origin in GUT scale physics. 
In this paper, we determine the values of GUT scale parameters which lead to models of 
natural SUSY. We find that, at the GUT scale, 

• third generation mass parameters, m.o(3) ~ 0.5 — 4 TeV, 

• first/second generation mass parameters, 7ti-o(1,2) ~ 5 — 25 TeV, 

• unified gaugino mass parameters, ini/2 ~ 0.3 — 1.7 TeV, and 

• the trilinear (third generation) scalar coupling, Ao/"^o(3) ~ —2 

yield models with a natural SUSY spectrum. The range of tan/3 and mA are relatively 
unrestricted. Note that there is an upper bound on mo(l) 2): values much larger than about 
25 TeV push third generation squarks into the tachyonic range via 2-loop RGE effects. We 
also find that values of mh ~ 125 GeV are very difficult to reconcile with a spectrum with 
very light third generation scalars (mg(3) < 0.5 TeV), but values of mh up to 124 GeV 
can be realized if we allow rnq{3) up to 1 — 1.5 TeV instead. Since some third generation 
squarks and charginos are rather light in natural SUSY, the constraint from BF{b — t- sj) is 
rather strong, but models can be found which are consistent with the measured branching 
fraction. We provide some representative benchmark points for low and high values of mh. 
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Figure 11: Thermally averaged neutralino annihilation cross section times relative velocity in 
limit as w — > versus r?zo(3) for NS models with /i = 150 GeV together with the bound from the 
Fermi satellite on the cross section times velocity for WIMP annihilation to M^-pairs (assuming 
higgsino-like WIMPs saturate the measured dark matter abundance). 



At the LHC, the higgsino-like light charginos and neutralinos have only small energy 
release in their decays, and so will be difficult to observe, as noted for the related "hidden 



SUSY" scenario [47 1 . However, in the case of natural SUSY, all four third generation 
squarks may be produced at observable rates. Sometimes, the lightest of these may have 
just one decay mode accessible {e.g. case NSl in this paper), and thus may be described 
by an analysis using simplified models (this is essentially impossible if 5i is the lightest 
squark). However, other cases arise where several different cascade decay possibilities are 
open. The heavier third generation squarks decay via numerous modes, and could lead 
to novel signatures involving h or Z from their cascade decays. Gluinos are favored to be 
rather heavy and frequently beyond LHC reach, although cases where rrig ~ 1.5 TeV do 
occur especially for mq;(3) ~ 1 — 1.5 TeV. 

At a linear e^e~ collider, we expect pair production of the higgsino-like light charginos 
and neutralinos to offer a lucrative discovery program of physics though specialized search 
strategies will be needed to pull out the rather soft signal events. In addition, it is possible 
that several third generation squarks and sleptons may be accessible to a LC with y/s 
extending up to ~ 1 TeV or beyond. Although these may decay to the light chargino as 
well as two lighter neutralinos, it will be challenging to sort out the various signals from 
the electroweak-ino cascades with very small secondary mass gaps. To our knowledge there 
are no dedicated studies for event topologies with this novel spectrum. 

In Natural SUSY models, the higgsino-like neutralino Zi is lightest MSSM particle. 
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and standard relic density calculations predict an underabundance of higgsino-like WIMPS 
by a factor typically 15. Such an under-abundance can be easily boosted to higher values if 

1. there are late decaying moduli fields with large branching fractions to SUSY particles, or 

2. if the PQ solution to the strong CP problem is invoked, whereupon thermal production 
of heavy axinos followed by late-time decays in the early universe can augment the higgsino 
abundance. In this latter case, any remaining under-abundance can be filled by axions. In 
the case of higgsino dominance of the dark matter abundance, then we expect higgsino- 
like WIMPs to be detected by the next round of direct and indirect WIMP detection 
experiments. An axion signal could also be a viable possibility. 
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